Thermal switching rate of a ferromagnetic material with uniaxial anisotropy 
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The field dependence of the thermal switching rate of a ferromagnetic material with uniaxial 
anisotropy was studied by solving the Fokker-Planck equation. We derived the analytical expression 
of the thermal switching rate using the mean first-passage time approach, and found that Brown's 
formula [Phys. Rev. 130, 1677 (1963)] is applicable even in the low barrier limit by replacing the 
attempt frequency with the proper factor which is expressed by the error function. 
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I. INTRODUCTION 

The thermal switching rate is an important quantity 
of ferromagnetic materials for applications in spintron- 
ics devices such as a magnetic recording media and spin 
random access memory (Spin RAM). The smaller the 
thermal switching rate, the longer the data retention pe- 
riod. In most experiments^— , the thermal switching rate 
is obtained by measuring the magnetization switching 
probability in the thermally activated region with the 
assistance of the applied field H app \ or electric current / 
(spin torque^) . Then the thermal switching rate can be 
analyzed by Brown's formula^ which is given by 
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where a, 7, and 7?k are the Gilbert damping constant, 
gyromagnetic ratio, and uniaxial anisotropy field, re- 
spectively Ao = M HkV I '{2k-QT) is the thermal sta- 
bility, where M , V , and T are the magnetization, vol- 
ume of the free layer, and temperature, respectively. 
h = —H app i/HK is the ratio of the applied field to the 
uniaxial anisotropy field, and Ao(l — h) 2 is the barrier 
height of the magnetic free energy. 

Brown's formula was derived in the high barrier limit 
Ao(l — h) 2 S> 1. In this case, it takes a long time to 
observe the thermal switching (more than a day or a 
year, depending on the value of Ao). On the other hand, 
the most interesting and important limit in the experi- 
ments is the low barrier limit with high thermal stability 7 
because the experimental determination of the switching 
rate, as well as the thermal stability, requires a large num- 
ber of observations of switching: For example, in Rcf. 2 , 
4000 times of switching were observed for one sample. 
Thus, in the experiments, a high field \h\ < 1 is ap- 
plied to the ferromagnetic materials to quickly observe 
the thermal switching. In such limits, Brown's formula 
gives unphysical predictions, as shown in this paper, and 
is no longer applicable. However, Brown's formula has 
been widely used in experiments to determine the switch- 
ing rate and thermal stability^—, which may, for exam- 
ple, lead to a significant error in the evaluation of the 
retention period of Spin RAM. Thus, it is important to 
study the validity of Brown's formula in the low barrier 



limit with high thermal stability. Also, the derivation of 
a simple and useful analytical expression of the thermal 
switching rate is desirable for the experiments. 

In this paper, we studied the field dependence of the 
thermal switching rate of ferromagnetic materials with 
uniaxial anisotropy by solving the Fokker-Planck equa- 
tion. We investigated the analytical expression of the 
mean first-passage time, and found that Brown's formula 
is applicable even in the low barrier limit by replacing 
the attempt frequency with the proper factor, which is 
expressed by the error function. We also compared the 
analytical formulas in several limits with the numerically 
calculated values and confirmed the validity of each an- 
alytical formula. 

The paper is organized as follows. In Sec. [TTl we 
derive the theoretical expression of the switching rate 
of the uniaxially anisotropic ferromagnetic material by 
the mean first-passage time approach, and show that the 
mean first-passage time approach is reduced to Brown's 
formula in the high barrier limit. Sections Mil and IIVI are 
the main parts of this paper. We derive the analytical 
formulas of the switching rate in the low barrier and low 
thermal stability limits and compare these to the numer- 
ically calculated switching rate. Section [V] is devoted to 
the conclusions. 



II. MEAN FIRST-PASSAGE TIME OF 
MAGNETIZATION SWITCHING 

We first derive an analytical expression of the mean 
first-passage time of the ferromagnetic material with uni- 
axial anisotropy by solving the Fokker-Planck equation. 
We assume that the dynamics of the magnetization in the 
ferromagnetic material is aptly described by the Landau- 
Lifshitz-Gilbert equation given by 
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= 7m x H 
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where m = (sin 9 cos tp, sin 9 sin ip, cos 9) is the unit vector 
pointing to the direction of the magnetization M. The 
magnetic field H = -dF/d(MV) = (H appl + H K cos6»)e z 
consists of the applied field and the uniaxial anisotropy 
field. F = -MH appl Vcos6 - (MH K V/2) cos 2 9 is the 
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magnetic free energy. In the uniaxially anisotropic sys- 
tem, the magnetization dynamics is described by 9 alone. 
We assume that \h\ < 1 because we are interested in 
the thermally activated region. Then, the free energy 
has two local minima at 9 — 0,7r and a maximum at 
9 m = cos -1 h. The difference of the free energies at 
9 = and 9 = 9 m (i.e., the barrier height) is given by 
Ao(l — h) 2 . Since below we focus on the thermal switch- 
ing from 9 = to 9 = ir, we assume that the field is 
applied to the — z direction (i.e., h > 0). 

At finite temperature, the thermal fluctuation gives ad- 
ditional torque on the magnetization, which is described 
by —7m x h. Here the random field h satisfies the 
fluctuation-dissipation theorem, 
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where i, j = x,y, z and T is the temperature. (• • • ) means 
the statistical average. The stochastic motion of the mag- 
netization due to the thermal fluctuation is described 
by the probability function P(9, t\9' , t'), which represents 
the transition probability from the state 9' at time t' to 
the state 9 at time t. From Eq. (J2J, the Fokker-Planck 
equation for the probability function is given byS. 
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or in terms of W = Pj sin 9, 
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where 7' = 7/(1 + a 2 ). Brown 5 calculated the switching 
rate by approximately solving Eq. (|5|) or by investigating 
the smallest nonvanishing eigenvalues, which characterize 
the time relaxation of the probability function. In this 
paper, we employ the backward Fokker-Planck approach 
to derive the mean first-passage time approach, which is 
the same as the eigenvalues of Eq. (J5J) and is useful to 
evaluate the thermal switching rate in several limits. 

The mean first-passage time T(9')&&, which character- 
izes how long the magnetization in the potential F stays 
within the region < 9 < 9 m , is defined as 

T(9') = dt d9P(9,t\9 , ,0). (6) 
Jo Jo 

The equation to determine the mean first-passage time 
is obtained from the backward Fokker-Planck equation 10 
as 
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and is given by 
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To solve Eq. (JSJ), it is convenient to transfer the variable 
6' to z = cos 9'. Then, Eq. ((8]) can be expressed as 
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The mean first-passage time is obtained by solving Eq. 

with the appropriate boundary conditions. We use 
the reflecting and absorbing boundary conditions at 9' = 
and 9' = 9 m , respectively^, that is, (dT/d9') g / =0 = 
and T(0 m ) = 0. Then the mean first-passage time is 
given by 
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Here, x and y are the integration variables. Equation (|I0[) 
was derived in RefJ^. Once the magnetization reaches 
9 = 6> m , it moves to 9 < 9 m or 9 > 9 m with the probability 
1/2. Thus, the switching rate is given by&^i 
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Below, we denote T(z = 1) as T. For the later discussion, 
we also define the normalized (dimensionless) switching 
rate as 



2A 



■i?. 



It should be noted that the y integral in Eq. (fTUf can be 
written in a different form by introducing the imaginary 
error function erfi(z) = (2/^/tt) J* dye y , that is, 

dyexp [A {y- h) 2 } 

crfi (v / Ao(l - hj\ - crfi Cy/K^(x - h)^j . 

(13) 

By using the Taylor expansion of the imaginary error 
function around x = 1, we found that the mean first- 
passage time can be expressed as 
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where the coefficient Ci = C^(Ao, h) is given by 
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The first three terms of Ce( Ag, h) are given by 

Ci=e A (l-k)^ (16) 

C 2 = -A (l-/i)e Ao ( 1 - /l ) 2 ! (17) 

2] „A (l-h) 2 _ Qg) 



C 3 = ^[l + 2A (l-/ l ) 2 ]e 

References&ii~— also calculate the switching rate by 
analytically or numerically solving the forward or 
backward Fokker-Planck equation. For example, in 
Refs.— i 12 i 15 i 18 i 19 ; the switching rate is calculated by ex- 
panding the probability function with the Legendre poly- 
nomials. In RefAi, the mean first-passage time after 
the integration shown in Eq. (|13l) is derived in terms 
of the Kummer's function, M(a, b, z) = Y]^°_ \T(a + 
n)T{b)z n ]/[nW(a)T(b + n)], where T(z) = J °° dtt^e^ 
is the r function [the Kummer's function satisfies the re- 
lation erfi(tt) = (2w/ v /7r)A/(l/2,3/2,it 2 )]. To obtain the 
analytical formula of the switching rate, these studies 
mainly investigated the high barrier (Ao(l — h) 2 ^> 1), 
small applied field (h ~ 0), or low thermal stability 
(Ao <C 1) limits. Although these constitute exact so- 
lutions, the particular limit of high fields and high ther- 
mal stability treated in the present manuscript was not 
singled out for particular attention, apart from a few nu- 
merical evaluations of the mean first-passage time in this 
regimeii. On the other hand, we note that the expansion 
shown in Eq. (|14[) is useful in deriving the analytical for- 
mula of the switching rate in the limits of Aq(1 — h) 2 <C 1, 
h < 1, and Ao 3> 1, which can be directly applied to an- 
alyze the experiments to evaluate the thermal switching 
rate and thermal stability accurately. 

Before proceeding with further calculations, let us give 
a brief comment on the spin transfer torque. In a fer- 
romagnetic multilayer with a pinned layer, the electric 
current applied to the multilayer exerts a spin trans- 
fer torque^, which gives the additional term jH s m x 
(p x m) to Eq. ©. Here H s = fir]I/{2eMV) and 
p are the strength of the spin transfer torque and the 
unit vector pointing in the direction of the magnetiza- 
tion of the pinned layer, respectively. The positive cur- 
rent (/ > 0) with the spin polarization rj is defined as 
the electron flow from the pinned to the free layer. In 
general, the effect of the spin torque cannot be included 
into the magnetic free energy, and the steady-state so- 
lution of the Fokker-Planck equation deviates from the 
Boltzmann distribution. Thus, the theoretical approach 
to the thermal switching rate shown in Ref. 6 is no longer 
applicable. However, in the uniaxially anisotropic system 
with p |j e z , the effect of the spin transfer torque can be 
taken into account by replacing H app \ in the magnetic 
free energy by i? apP i + H s /o£t—. Thus, the following 
discussions are applicable to the spin transfer torque sys- 
tem in this special case. It should be noted that Suzuki 
et al£^ showed that the results in the special case can be 
applied to general systems by replacing the parameters 
(for example, H^) with the appropriate values. 

At the end of this section, let us show that Eq. (fTU)) 



reproduces Brown's formula in the high barrier limit, 
A (l — h) 2 » 1. In this limit, the integral in Eq. 
(|10[) is dominated by the contribution around x = h. 
Since lim 2 _>o erfi(z) = 0, the result of the y integral [Eq. 
(fT3|)] in the limit of x — > h is given by v / 7rerfi[ v / Ao(l — 
h)]/(2y / Ao). By using the following formula, 



dye y 



we found that 
1 
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On the other hand, the x integral is given by 

1 e -A {x-h) 2 1 rl-h 

dx 
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where erf(z) = (2/i/tF) J* dye~ y is the error function. 
Since lim.2_j.00 erf (2) = 1 in the high barrier limit, the 
x integral can be approximated to y / 7r/[2 v / Ao(l — h 2 )]. 
Then, Eq. ([XDJ yields 
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The switching rate i?i = 1/(271) is identical to Brown's 
formula, Eq. ((T|). We can easily find that R\ is zero in 
the high field limit (h — > 1). Also, the attempt frequency, 
which is defined by Rie Ao ^~ h:> , increases with decreas- 
ing temperature in the limit of T — > 0. However, physi- 
cally, the switching rate should increase with increasing 
field strength, and the attempt frequency should be zero 
in the zero temperature limit. The origin of these contra- 
dictions is the high barrier assumption, Ao(l — h) 2 3> 1, 
in its derivation. 



III. ANALYTICAL EXPRESSIONS OF 
SWITCHING RATE IN OTHER LIMITS 

In this section, we derive the analytical expressions of 
the switching rate in the low barrier [Ao(l — h) 2 <C 1] 
and low thermal stability (Ao < 1) limits. 

First, we investigate the analytical expression of the 
switching rate in the high field limit (h —> 1) with a high 
thermal stability (Ao 3> 1), that is, the low barrier limit 
[Ao(l — h) 2 <C 1]. In this limit, it is sufficient to take 
into account the term of £ = 1 in Eq. (fT_j) , and the mean 
first-passage time is given by 



T 2 
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Although Ci for I > 2 is larger than unity for Ao ^> 1, 
we can neglect the higher order terms of Ci in Eq. ([14"]) 
because the integral interval, [h, 1], is very close to x — 1, 
and thus, the factor (1 - x) 1-1 for I > 2 in Eq. ([13} can 
be approximated to zero. Similarly, the factor 1/(1 + x) 
in the x integral can be approximated to 1/2. Then, the 
switching rate i?2 = 1/(272) is obtained as 



R2 



-Ao(l-ft) 2 



erf [VA^(1 - h)} ' 



(24) 



This is the analytical expression of the switching rate 
in the low barrier (high field) limit, and is one of the 
main findings in this paper. It should be noted that i?2 
increases with increasing field, and the attempt frequency 
R 2 e Aotyl ~ h ^ vanishes in the limit of T — > 0, which agrees 
well with our intuition. Equation ([24]) can be applied to 
the experiments to evaluate the thermal switching rate 
i?2 and the thermal stability Ao. 

We also derive the analytical expression of the switch- 
ing rate in the low thermal stability limit (A <C 1). 
In this limit, Cg in Eq. (fl"5"]) can be approximated to 
Ce ~ 8i t x. Then, Eq. (fl4l can be approximated to 



r 3 



2A 



dx- 
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Thus, the switching rate R% = 1/(273) is given by 



Rs 



4A log[2/(l + h)Y 



(25) 



(26) 



The switching rate R3 increases with increasing field 
magnitude, and vanishes in the limit of T — > 0. It should 
be noted that Eq. ([2"6"]) is applicable to the whole range 
of the field, \h\ < 1. We also noted that Eq. ([26]) is 
consistent with the result of Klein derived using a differ- 
ent approach^. Coffey et 



M- also derived Eq. ([26]) in 
which the concept of the mean first-passage time for a 
spherical domain was first introduced. 

By comparing Eqs. ([!]) and ([24]). we found that the 
switching rate is given by the product of the attempt fre- 
quency and the exponential factor e - A o( 1 -' 1 ) j anc l that 
Brown's formula is applicable even in the low barrier limit 
by replacing the attempt frequency, a.'y'H^l — h 2 ){l — 

/A^(l- 



h)y/Ao/%, by cry'H K /[y/TrA eT{(y/^(l-h))]. However, 
such exponential dependence is valid only in the high 
thermal stability limit: In the low thermal stability limit, 
the field dependence of the switching rate is described by 
the logarithm function log(l + h), as shown in Eq. ([2"6")> . 



IV. COMPARISON WITH NUMERICAL 
CALCULATION 

Finally, we compare the analytical formulas of the 
switching rate [Eqs. ([I]), flU}, and (|2"6)] with the nu- 
merically calculated values of Eq. (ITU1) . For comparison, 
it is convenient to evaluate the normalized switching rate, 
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FIG. 1: The field h dependences of the normalized switching 
rates, f and rk, where r is the numerical value of Eq. (|12|) 
while n, ri, and rz are Brown's formula [Eq. I|27[l]. high field 
limit [Eq. (T28l) ]. and low thermal stability limit [Eq. ((29)) ]. 
respectively, (a), (b) f (solid), n (dashed), and T2 (dotted) 
in the low and high field region with Ao = 60, respectively, 
(c), (d) f (solid), n (dashed), and r$ (dashed-dotted) in the 
low and high field regions with Ao =0.1, respectively. 



(b) 




I -20 



0.2 0.4 0.6 0.8 1.0 
normalized applied field, h 



0.2 0.4 0.6 0.8 1.0 
normalized applied field, h 



FIG. 2: The field h dependences of the errors of the ana- 
lytical formula, et- (a) ei (dashed) and e2 (dotted) with the 
high thermal stability, (b) e\ (dashed) and e% (short dashed- 
dotted) with the low thermal stability. 



Eq. lfl"2"]) . Due to the normalization, the switching rates 
depend on Ao and h only. We denote the numerical value 
of Eq. ([T2"]) as r, and the normalized values of Eqs. ([1}, 
([M]) . and ([2"6"]) as r\, r 2 , and r 3 , respectively. Explicitly, 
these are given by 



ri 



'=± (l-h 2 )(l-h)e- A ^- h ^, (27) 
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The normalized values of r-i and r 3 in the limit of h — > 1 
are identical, 1/(1 — h). For the quantitative discussion, 
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we also define the error of the analytical result by e k = 
[{r k - r)/f] x 100 %. 

Figures Q] (a) and (b) are the field dependences of r 
(solid), r\ (dashed), and r 2 (dotted) with the high ther- 
mal stability Ao = 60 in the low and high field regions, 
respectively. As shown, n and r 2 almost overlap with the 
exact result (f) in the low and high field regions, respec- 
tively. For example, the errors of Brown's formula (r\) 
are e x = 17.7 and -88.4 % for h = 0.5 and 0.95 while 
those of r 2 are e 2 = —94.8 and —4.5 % for h = 0.5 and 
0.95, respectively: see Fig. H](a). The good agreement of 
our formula [Eq. (|24[) ] in the high field limit is important 
for the experiments to guarantee the validity of Brown's 
formula. Moreover, Brown's formula gives unphysical re- 
sults in the limit of h . — > 1; that is, the switching rate 
decreases with increasing field, as shown in Fig. Q] (b), 
while Eq. (|24[) gives a physically reasonable and quanti- 
tatively good estimation of the switching rate. 

In Figs. [T](c) an d (d), we show the normalized field de- 
pendences of r (solid), r\ (dashed), and rj, (short dashed- 
dotted) with the low thermal stability Ao = 0.1 in the 
low and high field regions, respectively. The errors, e\ 
and e3, are shown in Fig. [2] (b). One can clearly see that 
r3 agrees well with the numerical result: e% = 0.7 and 
0.007 % for h = 0.5 and 0.95, respectively. On the other 
hand, Brown's formula reduces e\ = —99.2 and —99.9 % 
for h — 0.5 and 0.95, respectively. The large difference 
between f and r\ is due to the break down of the high 
barrier assumption in Eq. (fl}. 

Let us briefly discuss the application of our formula 
[Eq. (1241) ] to spintronics applications. Since the switch- 



ing rate estimated by Brown's formula is much smaller 
than the exact value, as shown in Fig. [1] (b), a relatively 
small value of the thermal stability as a fitting parame- 
ter is required to fit the experiments with Eq. (fTJ) . Then, 
the estimated value of the retention time of Spin RAM, 
which is proportional to exp(Ao), is significantly under- 
estimated. On the other hand, Eq. (|2"4")l can give accu- 
rate estimates of the thermal switching rate, the thermal 
stability, and also the retention time of Spin RAM. 



V. CONCLUSIONS 

In conclusion, we studied the field dependence of the 
thermal switching rate of a ferromagnet with uniaxial 
anisotropy theoretically. We derived the analytical ex- 
pression of the switching rate by the mean first-passage 
time approach, and showed that Brown's formula is ap- 
plicable even in the low barrier limit by replacing the 
attempt frequency with the proper factor which is ex- 
pressed by the error function. We also compared the 
analytical formulas of the switching rate for several lim- 
its with the numerically obtained values and showed the 
validity of each analytical formula. 
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